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ABSTRACT 



Context. The space missions CoRoT and Kepler provide high quality data that allow us to test the transport of angular momentum in 
stars by the seismic determination of the internal rotation profile. 

Aims. Our aim is to test the validity of the seismic diagnostics for red giant rotation that are based on a perturbative method and to 
investigate the oscillation spectra when the validity does not hold. 

Methods. We use a non-perturbative approach implemented in the ACOR code (Ouazzani et al. 2012) that accounts for the effect of 
rotation on pulsations, and solves the pulsation eigenproblem directly for dipolar oscillation modes. 

Results. We find that the limit of the perturbation to first order can be expressed in terms of the rotational splitting compared to the 
frequency separation between consecutive dipolar modes. Above this limit, non-perturbative computations are necessary but only one 
term in the spectral expansion of modes is sufficient as long as the core rotation rate remains significantly smaller than the pulsation 
frequencies. Each family of modes with different azimuthal symmetry, m, has to be considered separately. In particular, in case of 
rapid core rotation, the density of the spectrum differs significantly from one m-family of modes to another, so that the differences 
between the period spacings associated with each ra-family can constitute a promising guideline toward a proper seismic diagnostic 
for rotation. 

Key words, asteroseismology - stars: interiors - stars: oscillations 



1. Introduction 

Seismic measurements of rotation profiles inside the Sun as well 
as stars provide tight constraints on models of transport of angu- 
lar momentum (| Pinsonneault et al. 1989t lZahnll 992: Zah n et alJ 
' ll997tlTalon & Charbon nell2008l and references therein). In par- 
ticular, stars in late stages of evolution, due to the highly con- 
densed core, oscillate with non-radial modes that have a mixed 
character: they behave as p modes in the envelope and as g 
modes in the core. These modes, also known as mixed modes, 
are of particular interest for the determination of the rotation 
profile throughout the star, as they carry the signature of the 
star's inn ermost layers and are detectable at the surface. The 
CoRoT dBaglin et alJ l2006h and Kepler dBorucki et al.l l2010h 
spacecrafts have dramatically improved the quality of the avail- 
able asteroseismic data. Several recent studies reported the de- 
tection of mixed mode s that are split by rotation in a subgiant 
dDeheuvels et al.l I2012I) and in several red giants (Bec k et all 
2012; Mo sser et al.l l2012blal) observed with Kepler. A large 
number of these stars exhibit frequency spectra that show a quite 
simple structure where symmetric patterns around axisymmet- 
ric modes are easily identified. They are interpreted as multi- 
plets of modes split by rotation. The rotational splittings, i.e. 
the frequency spacing related to the lift of degeneracy caused 
by rotation is then used to determine the core rotation. The val- 
ues of the corresponding splittings are quite small and the use 
of the lowest order approximation to derive the splittings from 
stellar models can be justified. Such studies led to the determi- 
nation of unexpectedly low central rotation frequencies (of few 



hundreds of nHz). These results are in strong disagreement with 
the core rotation frequencies predicted by ev olutionary models, 
which are of the order of fe w dozens of //Hz dEggenberger et all 
120121 ; iMarques et al.l 120131) . They show that the transport pro- 
cesses currently included in stellar models are not able to spin 
down the core of red giant stars enough to explain the slowly 
core rotating red giants. On the other hand , a large set of red gi - 
ant stars show complex frequency spectra (Mosser et al. 2 012al) . 
in particular with non symmetric multiplets and therefore are 
likely rotating fast. Their rotation must then be investigated with 
non-perturbative methods. 

In this context, we first report on the relevance of using a 
first order approach for the inference of rotation from seismic 
spectra of red giant stars with slowly to rapidly rotating core. 
When not relevant, we adopt the non-perturbative approach in 
order to shed light on the behaviour of splitted mixed modes in 
red giants spectra. 



2. Theoretical frequency spectra for red giants 

This work is based on the study of a model, (model Mi with 
mass 1.3 M and radius 3.6 R ), at the bottom of the red gi- 
ant branch. The stellar model is computed with the CESTAM 
code (Code d' Evolution Stellair e, avec Transport, Adaptatif et 
Modulaire, Marques et al. 2013). Transport of an gular momen - 
tum induced by rotation is included according to IZahnl (Il992l) . 
The central rotation rate in this model is Q c /27r ss 180//Hz, while 
the surface convective region rotates at a rate of 1 yuHz. The 
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Fig. 1. Pulsation frequencies in the inertial frame v (in yuHz) 
versus core rotation frequency Q. c /2n (in pMz) for dipolar mul- 
tiplets. Red curves refer to m = -1 (prograde) modes, green 
curves to m = modes (axis-symmetric) , and blue curves to 
m = 1 (retrograde) modes, for central rotation frequency rang- 
ing from to 100 fiHz. 
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Fig. 2. Top: Ratio between the contributions to kinetic energy of 
the i — 3 component and the I = 1 component of a dominated 
£ = 1 mode. Bottom: Difference between the pulsation frequen- 
cies computed with three coupling terms (M = 3) and with one 
coupling term (M — 1), divided by the frequency separation be- 
tween two consecutive mixed modes. Computed for the triplet 
indicated in thick lines in Fig. 



distortion of the model due to the centrifugal force can be ne- 
glected e verywhere in red giant star s. For the model Mi, scal- 
ing laws (Kiel dsen & Bedding]! 19951) give a frequency of max- 
imum power around v max = 289/zHz and a large separation 
of Av = 23 juHz. We then compute frequencies ranging be- 
tween v max + 2Av. In this frequency range, the impact of the 
Coriolis force remains small except in the very inner layers of 
the star where it can significantly affect the modes. In order to 
investigate the effect of core rotation on the frequency spec- 
trum, we compute sets of frequencies for model Mi for a se- 
quence of rotation profiles. This sequence is obtained by divid- 
ing the whole rotation profile given by CESTAM for model Mi 
by constant factors. The oscillation frequencies are calculated 
by the non-perturbative pulsation co de ACOR (Adiabatic Code 
of Oscillations including Rotation, lOuazzani et al]|2012l) . The 
eigenmodes are obtained as a result of the coupling of spherical 
harmonics. In what follows, M coupling terms means expansion 
on the t = 1,3,..., 2M+ 1 spherical harmonics for the scalars and 
the poloidal velocity component and i = 2,4,..., 2M + 2 for the 
toroidal velocity component. 

Figure Q] shows the frequencies of several dipolar (I - 1) 
multiplets with increasing rotation frequency, computed using 
one coupling term (M = 1) in the spectral expansion. Starting 
at low rotation rate, due to the combined action of the Doppler 
effect and the Coriolis force, in the inertial frame, prograde 
(m — -1) modes are shifted towards higher frequencies, whereas 
retrograde (m = +1) modes are shifted towards lower frequen- 
cies. Crossings between modes of different symmetry occur ap- 
proximately at £l c /2n > 8//Hz (for model Mi). As a result, for 
core rotation frequencies above 8yt/Hz, modes of different m are 
no longer gathered by original multiplets, i.e. triplets of modes 
which have the same degenerate frequency without rotation. 

The choice of M — 1 coupling term in the spectral expansion 
is indeed sufficient for most of the rotation profiles investigated 
here. This is illustrated in Fig. |2]which shows the impact of the 
number of spherical harmonics used in the eigenmodes expan- 
sion. Fig. [2] present the results of non-perturbative calculations 



using three coupling terms (M = 3). Most modes are clearly 
dominated by their t = 1 component (Fig.|2l top). Figure ^(bot- 
tom) shows that the frequency difference between calculations 
including one and three spherical harmonics are smaller then the 
frequency difference between two consecutive mixed modes (for 
Q. c /2n < 100;uHz by two orders of magnitude). From now on, 
we adopt M — 1 . 

The trapping of mo des essentially de pends on their fre- 
quency (see for instance lUnno et al. 1989). As seen in Fig. Q] 
when Q. c increases, so does the frequency difference between 
the three members of a triplet and the trapping of these mem- 
bers can be significantly different, to such an extent that they 
end up with very different p-g nature. This is illustrated in Fig. 
[3] where the kinetic energy corresponding to modes of different 
m around v max is plotted. It gives an indication of the p-g na- 
ture of modes: the p-dominated modes (referred as p-m modes) 
correspond to the minima of energy, while the g-dominated ones 
(g-m modes) are associated with the maxima. The three modes 
circled in black belong to the same original multiplet, and show 
different p-g nature. The prograde and the retrograde modes are 
g-m modes, while the m = one is a p-m mode. This change 
of nature induced by rotation, which depends on the azimuthal 
order, can also be characterized by the number of nodes in the 
p-mode cavity (n p ) and in the g-mode ca vity (n g ) (calcu lated 
according to the Cowling approximation. ICowhngl 09411) . As 
shown in Fig. [3] the n p and n g values of the members of an orig- 
inal triplet are modified differently. These gains (for m — 0, 1 
modes) or losses (for m = -1 modes) of nodes occur during 
gravito-acoustic avoided crossings. When increasing the rota- 
tion rate, if the frequencies of two modes of same symmetry 
become very close, they avoid to cross each other and exchange 
nature. These avoided crossings need to be taken into account 
all the more because they mainly affect the p-m modes that are 
the most likely to be observed. Therefore, modes of different az- 
imuthal order m probe differently the stellar interior. Even if they 
belong to the same triplet with radial order n and degree (, they 
can be of very different nature, and using the rotational splitting 
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Fig. 3. For a central rotation frequency of 20//Hz, kinetic en- 
ergy of modes around a p-m mode as a function of the mode 
frequency. The three modes circled in black belong to the same 
original triplet. Inside the parentheses are indicated their number 
of nodes in the p cavity and in the g cavity. 

5u) n j = (<jj n ,t,m - w„/- m )/2m to determine the rotation rate is 
therefore questionable for red giants. 

3. Slow to moderate core rotating red giants 

For models with Q. c /2jt < 8 //Hz (see Fig. [TJ, the frequencies 
behave linearly with respect to Q c . In this range, the sectorial 
modes (f = | m |) are symmetrically distributed around the ax- 
isymmetric modes, and have the same radial order. In this case, 
a first order perturbation approach approximates well the ef- 
fects of rotation on the mode frequencies. The pulsation co n .i. m 
is then given by ui n ,t,m = <^n,t,o + rn5to l n ( ([Cowling & Ne wing 
1949 and lLedo"uxl [l949) where the first order rotational split- 
ting 6co l n . is expressed as a weighted measure of the star rota- 

tion rate: 8to n ( = L K„/(r) Q,(r) dr, where K n /(r) are the ro- 
tational kernels of the modes. They depend on the equilibrium 
structure and on the e igenfunctions of the unperturbed modes 
(see e.g. lGoupilll2009L and references therein). This formulation 
relies on two main assumptions. First, only the Coriolis force is 
significant and it is enough to account for its contribution to or- 
der Q(Q./ u)„.t). Second, it is derived fro m a variational principle 
(see e.g. iLvnden-Bell & Ostrikerl u967). which requires that the 
eigenfunction of the mode perturbed by rotation is close to the 
eigenfunction of the corresponding unperturbed mode. This is 
the case for model Mi for Q. c /2tt < 8 //Hz. 

Between 8 and 20 /iHz (Fig. [TJ, the mode frequencies seem 
to behave linearly with Q e , but are no longer gathered by orig- 
inal multiplet. Figure [4] (top) displays the apparent rotational 
splittings, taking half the difference between the closest modes 
of opposite azimuthal order m = +1. The values taken by the 
apparent rotational splittings range from few hundreds nHz to 2 
//Hz. This is of the order of the splittings measured in observed 
spectra. Because the apparent multiplets occur by 'accident', the 
curve in Fig. [4] (top) does not -and usually cannot- follow the 
particular 'V pattern around the p-m modes as ob served for 
some slowly rotating red g iants (see for observations Bec k et all 
l2012llMosser et ai1l2012bl) . This 'V pat tern has the same o rigin 
as the modulation seen for mode inerti a (|Goupil et aTJl2013l) that 
is due to the trapp ing of modes (see [Dziembowski et al. 2001; 
Dupret et al. 2009). The 'V pattern can then constitute a strong 
indication that the selection of multiplet has been done correctly 
in observed spectra. Figure [3]confirms that modes which belong 
to the same original multiplet (circled in black) are found in dif- 
ferent apparent multiplets, and show quite different p-g nature. 
FigE](bottom) compares the original splittings -i.e. half the dif- 
ferences between frequencies of modes from the same original 
triplet with opposite m- with the splittings given by the first or- 
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Fig. 4. For the model Mi and a central rotation frequency of 
20//Hz, Top: spacing between two closest modes of opposite az- 
imuthal order m = +1, i.e. apparent rotational splitting. Bottom: 
spacing between prograde and the retrograde modes which be- 
long to the same original multiplet (grey, crosses) and splittings 
computed with the first order perturbative method (green, plus 
signs). 



der approach. The g-m modes (e.g. at 277 //Hz and 303 //Hz) 
have the largest splittings. For such modes, both computations 
give very close results (with a discrepancy of ~ 4%). The p-m 
modes (e.g. at 290 //Hz and 315 //Hz) have the smallest splittings 
and show the expected 'V shape variation. For these modes, the 
discrepancies can reach 30%. This is due to the trapping of these 
modes which differs from one m to the other. Note that original 
splittings including three terms (M = 3) in the non-perturbative 
calculations give relative differences lower than 10~ 6 when com- 
pared with results of computations using two Y? (Fig. |2j. 

Therefore, for a rotation between 8 and 20 //Hz, the fact that 
the multiplets are correctly selected can be assumed if the split- 
ting given by apparent triplets do follow the 'V pattern. If so, 
the first order approach provides an order of magnitude esti- 
mate of the core rotation rate, but non-perturbative modeling is 
required for accurate quantitative conclusions on the rotational 
profile from p-m modes. 



4. Rapid core rotating red giants 

For Q > 20/zHz, the mode frequencies no longer behave lin- 
early with the rotation rate (see Fig. QJ. This is due not only to 
higher order effects that come into play, but also to the trapping 
of modes that is modified by rotation. Modes which belong to 
the same original multiplet do not have the same radial distribu- 
tion. Under these circumstances, the notion of rotational splitting 
as defined by 6oj l n , is no longer relevant, and cannot be simply 
related to the rotation profile. Figure |5]( top) shows the kinetic en- 
ergy of pulsating modes for a core rotation rate of 140 //Hz. The 
large separation, as measured by the difference between consec- 
utive p-m modes with minimum kinetic energies, is conserved. 
The higher kinetic energy of prograde modes indicates that they 
are more of g nature than retrograde or axisymmetric ones. This 
is due to the shift in frequency induced by rotation that brings, in 
the same frequency range, modes which belong to very different 
parts of the zero rotation spectrum (see Fig. [TJ. 

Figure [5] (top) shows that, in the same range of frequency, 
there are more m — -1 modes than m = modes, and in turn, 
there are more m — modes than m = +1 ones. In order to 
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Fig. 5. For a central rotation frequency of 140 //Hz, top: pul- 
sation modes' kinetic energy with respect to their frequency. 
Bottom: period spacing with respect to the pulsation frequency. 

highlight this difference in distribution, the period spacings are 
plotted for each to value in Fig. [5] (bottom). In this diagram, the 
families of modes of different azimuthal orders clearly show dif- 
ferent values of period spacing, ranging from the lowest for the 
prograde modes to the highest for the retrograde ones. Note that 
this phenomenon appears only for a rotation high enough that the 
distributions of the three TO-families of modes become clearly 
different (Q c ~ lOCty/Hz for model Mi). If one is able to mea- 
sure three different values of period spacings in a observed spec- 
trum, one is able to identify the values of the azimuthal order. 
Moreover, based on our study of synthetic spectra, we find that 
the differences between the values of period spacings increase 
with Q c , and therefore can constitute a guideline toward a con- 
straint on the rotation profile. 

5. Discussions and conclusions 

Some observed spectra of red giant stars show structures with 
nearly symmetric spacings around axisymmetric modes. For this 
kind of spectra, we found that one can figure out if these appar- 
ent multiplets correspond to the original splittings by plotting 
the apparent splittings as a function of the axisymmetric modes 
frequencies. If the resulting curve follows the "V" pattern, then it 
is reasonnable to assume that the triplets are correctly selected, 
and quantitative information on rotation can be derived from the 
observed splittings. This corresponds to a low rotation regime 
when the splitting is smaller than the frequency spacing between 
two consecutive modes of same to, i.e. for Q c /(27r)/(AP/P 2 ) < 2. 

If the apparent splittings do not follow the 'V pattern, the 
triplets overlap. The star is therefore located in the regime of 
moderate rotation. The frequency differences between the origi- 
nal triplet members are large enough that the trapping properties 
differ from one member to another. However, if it is still possible 
to sele ct correctly the original triplets (as done in Mosse r~et al.l 
1201 2al) . the first order approach could be used for a first guess of 
the core rotation rate. Non-perturbative calculations are never- 
theless required for a precise determination, particularly because 
of the avoided crossings that p-m modes undergo. This rotation 
regime corresponds to values f2 c /(27r)/(AP/P 2 ) between 2 and 5. 

Finally, some observed spectra of red giants do not show any 
regular or close to regular structures. We expect these spectra 
to correspond to rapid rotators. In this case, modes of different 
azimuthal orders have a very different nature. The concept of ro- 
tational splitting of modes with the same radial distribution and 



different to is no longer relevant, and one should consider sep- 
arately sub-spectra associated to each value of to. Provided that 
the rotation is high enough to give rise to clearly different dis- 
tributions with respect to to, the differences between the period 
spacings associated with to allows to identify the azimuthal or- 
der, and thereby offers the promising opportunity of deriving a 
proper seismic diagnostic. That corresponds to the very rapid ro- 
tation case i.e. for O c /(27r)/(AP/P 2 ) > 20. An intermediate case 
remains, where the rotational splitting is no longer a relevant 
seismic diagnostic, but where the rotation is too slow to allow 
to distinguish three different period spacings associated to the 
TO-families of modes (for Q c /(27r)/(AP/P 2 ) between 5 and 20). 
Establishing a diagnostic in this regime constitutes an important 
issue which needs to be adressed in a forthcoming paper. 

It can be surprising that when Q./(2jt)/(AP/P 2 ) < 20, on the 
one hand the perturbative approach gives inaccurate results, but 
on the other hand non-perturbative models where the spectral ex- 
pansions are limited to one term (plus the toroidal component) 
give accurate enough results. This can be understood with the 
help of the following equation: D(r)/27r ~ Av <k v max . In the 
frequency range of interest here, the rotation frequency remains 
significantly lower than v max throughout the whole star, and thus 
lower than the pulsation frequencies. Hence, the Coriolis terms 
are smaller than other terms in the equation of motion and the 
coupling of t = 1 modes with t = 3,5, ... components is small. 
However the rotational splitting is not small compared to the 
large separation Av. This implies that the trapping of each mem- 
ber of a same triplet can be very different, as illustrated in Figj3] 
Using the variational principle to model rotational splittings is 
thus not justified because the eigenfunctions are too different. 
The only way to properly model the effect of rotation on oscilla- 
tions is to solve separately the differential equations associated 
to each member of the multiplet. 

To sum up, the major impact of moderate rotation on red gi- 
ants' spectra is the modification of the trapping that depends on 
to. Let us precise that here the cavities are not considered to be 
modified by rotation, but it is the way modes probe them which 
differs from one member of a multiplet to another. In such a case, 
the first order perturbative approach gives inaccurate results and 
the TO-families of modes carry different information on the stel- 
lar interior. Only methods which compute these TO-sub-spectra 
independently, i.e. non-perturbative methods, are appropriate for 
studing the moderate to rapidly core rotating red giants. 

Acknowledgements. RMO thanks Benoit Mosser for fruitful discussions. The 
authors thank the referee for his comments that helped to improve the 
manuscript. 
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